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SUMMARY 

The primary objective of advanced STOL aircraft is the improvement of the 
nation's air transportation system by the elimination of delays and congestions 
associated with today's air travel. A new guidance technique, referred to as 
4-D guidance, is being developed to achieve this objective. The 4-D guidance 
technique synthesizes complex three-dimensional flight paths from a minimum 
set of input data and flies the aircraft along the paths according to a pre- 
specified time schedule. The two major elements of a 4-D guidance system are 
the trajectory synthesizer and the control law. Inputs to the trajectory synthe- 
sizer are the three-dimensional coordinates of way points, the turning radii, 
the speed ranges, the acceleration limits, and the arrival times at time control 
way points. First, the three-dimensional trajectory is computed by using circu- 
lar arcs and straight lines. Then the airspeed profile, compensated for wind, 
is calculated to achieve the desired arrival times. The synthesized trajectory 
is stored as a time sequence of reference states which the aircraft is forced to 
track by using a linear feedback law. 

INTRODUCTION 

If advanced STOL aircraft are to become viable elements of the nation's air 

transportation network, their unique performance characteristics must be exploited 

to the fullest possible extent. Most notable among these characteristics are the 
* 

aircraft's ability to make tight turns, steep climbouts and approaches, and to fly 
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over a wide range of airspeeds. Consequently, STOL aircraft will possess much 
greater flexibility in the terminal area than CTOL aircraft. On the other hand, 
the performance of STOL aircraft may be limited by two important factors: the 

aircraft's acceleration/deceleration capability, and the effects of wind. Some 
existing STOL aircraft are slow to decelerate at some speeds and thus require 
long deceleration periods in the terminal area. Since these aircraft will be 
able to fly at very low speeds, on the order of 60 to 65 knots, the effects of 
wind can become extremely important in predicting arrival times at various key 
points along the route. In light of these performance characteristics and limi- 
tations, the following questions have arisen: Can STOL aircraft follow very 

demanding ATC instructions? Can they fly complex three-dimensional trajectories 
in the terminal area? Are they able to maintain an accurate time schedule 
along their flight path? To answer these questions, a technique is needed to 
play the characteristics of any potential STOL aircraft against these terminal 
operating constraints. Such a technique will enable us to determine whether a 
potential STOL aircraft can operate within these constraints and, if not, to 
identify STOL aircraft deficiencies to aircraft designers. 

It seems likely that advanced short take-off and landing (STOL) aircraft 
will be equipped with guidance and navigation systems that satisfy far more 
stringent requirements than those currently used in conventional (CTOL) aircraft. 
Since these requirements are not yet clearly defined, investigations of the guid- 
ance problems for STOL aircraft under the most demanding circumstances were made. 
Although this approach places a significant burden on the aircraft and onboard 
guidance system, it ensures that the results obtained will not be invalidated by 
future developments in ground and airborne system design. 

It has been assumed that STOL aircraft will (1) fly curved climbouts and 
approaches using a microwave landing system; (2) operate in narrow airspace 
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corridors to avoid CTOL traffic, buildings, and other obstacles; and (3) per- 
form noise -abatement maneuvers. The guidance system that meets these require- 
ments is referred to as a 3-D guidance system, and it involves the precise 
description of complex curved flight paths and the control system to fly the 
aircraft along the paths. 

Precise path control alone, however, is not sufficient to achieve the 
principal goal of a STOL transportation system, which is the elimination of, or 
at least the substantial reduction in, undue delays and congestions. This goal 
can be accomplished by increasing the landing rates at airports; this increase, 
in turn, implies a greater precision in delivering aircraft to the runway thres- 
hold. In the current manual system the accuracy with which aircraft can be 
delivered to the runway is approximately ±15 sec. It has been suggested that 
this time can be reduced to about ±2 sec by spacing aircraft accurately with 
less vectoring and by minimizing pilot-controller communications. Both of these 
requirements can be fulfilled by achieving precise time control of the aircraft 
along its flight path. The guidance system that accomplishes this control is 
referred to as 4-D guidance. It is a guidance technique consisting of two major 
elements: the 4-D reference trajectory synthesis, and the control law to fly the 

aircraft along the reference trajectory. It must be emphasized that there is a 
very important, intrinsic difference between these two elements. By 4-D refer- 
ence trajectory synthesis is meant the generation of a precise 3-D flight path 
and a feasible speed profile along the path. Note that this process is a com- 
pletely open-loop, or predictive, process in the sense that it takes place before 
the aircraft arrives at the initial point of the trajectory. The action of the 
control law, on the other hand, is a closed-loop process going on in real time 



as the aircraft is tracking the reference trajectory. It is the combined open- 
loop— closed- loop process that is referred to as a 4-D guidance system. 

The purpose of this paper is to describe the 4-D guidance system that was 
developed at NASA Ames Research Center, and which will be flight tested in a 
STOL aircraft. Some of the crucial design considerations are discussed, and 
a brief description of the algorithms for trajectory synthesis and control law 
are given. Finally, some simulation results are presented to indicate the per- 
formance of the system. 

DESIGN CONSIDERATIONS AND INTERACTION WITH 
AIR TRAFFIC CONTROL 

Chief criteria for selecting the 3-D path are STOL terminal area maneuver ' 
requirements and simplicity of computation. These criteria are met by synthe- 
sizing the 3-D path from geometrically simple elements. In the horizontal 
plane these elements consist of segments of circles and straight lines. Com- 
plex flight paths are obtained by interconnecting several line segments and 
sections of circles with different radii. Paths constructed in this manner can 
yield minimum time trajectories as discussed in references 1 to 3. The vertical 
trajectory is synthesized from sections of constant flight-path angle. The 
complete three-dimensional flight path is then obtained by requiring the aircraft 
to fly the vertical profile along the ground track determined by the previously 
computed horizontal trajectory. A critical problem in implementing this syn- 
thesis procedure is minimizing the pilot work load in entering a trajectory. 

As explained in the next section, this problem is solved by using way points to 
specify the trajectory. 

After the three-dimensional path has been established, the desired position 
of the aircraft along the path as a function of time is determined from consid- 
erations of air traffic control. Generally, more than one aircraft will be 



flying along the 3-D flight path or will be merging with the path at certain of 
its points referred to as merging points. This effect is illustrated in 
figure 1. Aircraft (AC) on the two approach routes merge in the vicinity of the 
approach gate. In the current air traffic control (ATC) system the final approach 
controller is responsible for merging aircraft at this point. He does this by 
observation of the aircraft on the ATC radar and by issuing speed and vectoring 
instructions to the pilots. The delays, inaccuracies, and other limitations of 
this manual controller-pilot control loop yield a broad envelope of aircraft 
trajectories between the feeder fix and the gate, shown as the hatched area in 
figure 1. This manual technique can deliver aircraft to the approach gate with 
a time accuracy of approximately ±15 sec. It is well-known that the time control 
accuracy influences the required spacing on final approach, which, in turn, 
determines the landing capacity of the runway (ref. 4). Through the use of 4-D 
guidance techniques, the accuracy of time control can be greatly increased; thus, 
a basis for achieving higher landing rates and greater automation in the control ■ 
of terminal area traffic is provided. Instead of the controller issuing vector- 
ing and speed commands to space aircraft, the air traffic control system speci- 
fies only the desired arrival times at a small number of points along the termi- 
nal approach route and leaves the burden of computing aircraft control inputs to 
achieve these times to the onboard system. 

From a knowledge of the 3-D path and the desired arrival times at specified 
points on the path, the 4-D guidance system computes the required airspeed pro- 
file. The airspeed computation algorithm must consider the minimum and maximum 
permissible airspeed, the aircraft's acceleration and deceleration capability, 
the landing approach speed, and the effect of winds. Also, in specifying arrival 
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times, air traffic control will use only limited knowledge of the aircraft's 
performance capabilities. Therefore, the algorithm must first determine the 
feasibility of the specified arrival times by comparing them with the true 
minimum and maximum times the aircraft can achieve without deviating from the 
3-D path. If the specified arrival times cannot be achieved, air traffic con- 
trol must be requested to reassign arrival times or permit delaying maneuvers 
such as holding and/or path stretching. 

In the preceding discussion of arrival time assignment, the question left 
unanswered was how air traffic control will generate the arrival times to be 
assigned to each aircraft. To clarify this question, consider an aircraft 
equipped with a 4-D guidance system which has just arrived at one of the feeder 
fixes in figure 1. It is assumed that the aircraft had previously been cleared 
to proceed toward the feeder fix. Those aircraft currently flying between the 
feeder fixes and the runway were previously assigned arrival times at the gate 
and the touchdown point. Given the schedules of these aircraft and, for the 
new aircraft, an estimate of the minimum and maximum times to the gate and from 
the gate to the touchdown point, a technique is needed to specify feasible times 
at the gate and at the touchdown point such that separation standards between 
aircraft are satisfied and that the aircraft lands in minimum time. A general 
algorithm to calculate such times has been described in reference 5. 

SPECIFICATION OF 3-D PATH 

The problem of specifying and calculating the 3-D path is divided into two 
problems solved in sequence. First the projection of the trajectory in a hori- 
zontal plane is computed from an analysis of the way-point coordinates and the 
desired turning radii. Then the known arc length of the horizontal trajectory 
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together with the altitude difference between adjacent way points is used to 
determine the flight -path angle and, therefore, the altitude profile between 
adjacent way points. 

A crucial part in the calculation of horizontal trajectory parameters is 
the interpretation of way points. This part is explained with the help of an 
example trajectory shown in figure 2. The trajectory begins at lift-off and 
terminates at touchdown, way point 7 (WP7). Although this trajectory can be 
flown by a STOL aircraft, its shape has no significance beyond illustrating 
the construction procedure. All parts of the trajectory consist of segments 
of straight lines and circles. The basic problem to be solved can be stated 
as follows: What is the essential input information that the pilot must pro- 

vide to the airborne computer in order to generate this trajectory uniquely? 
The solution lies in the definition of two types of way points, which are 
referred to as the ordinary and final heading way points. 

Ordinary Way Points 

An ordinary way point is exemplified in figure 2 by way points, 2, 4, and 

5 (WP2, WP4, and WP5) . Its location is defined by the intersection of two 

straight lines. The two lines are connected by an arc of a circle tangent to 

both lines. Thus, the sharp corner at the way point is rounded to obtain a 

trajectory the aircraft can fly. The radius of the circle used in rounding 

the corner can either be explicitly specified by the pilot as in this example 

or it can be implicitly determined from a maximum bank angle constraint <J> . 

For a given , the minimum turning radius R . depends on the maximum 
6 max & mm r 

ground speed V that can be attained in a 360° turn and is 

6 r g,max 


R . 


mxn 


V 

g,max 
g tan ip 


max 


.CD 
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where g is the acceleration of gravity. The maximum ground speed is the sum 
of the maximum airspeed and the magnitude of the wind vector. 

Throughout this paper, the runway-centered coordinate system is used to 
specify points in the plane. The origin of this system is at the touchdown 
point, way point 7, the X-axis points in the landing direction and the Y-axis 
points to the right when facing in the landing direction. Heading angles are 
measured clockwise from the direction of the positive X-axis. 

Suppose the pilot has entered the x- and y-coordinates of a sequence of 
ordinary way points together with the turning radii to be used in rounding the 
corners at way points. From this information the onboard computer calculates 
various parameters defining the trajectory. It is convenient to compute these 
parameters successively, starting with the last way point to be flown through 
and ending with the first one. For this purpose figure 3 shows the trajectory 
between way points i - 1 and i + 1 and also defines various quantities used in 
the calculation. It is assumed that the calculations from way point i + 1 to 
the last way point have been completed. These calculations yielded the quanti- 
ties Z i+1 , x i+1 , x pi+1 , and i|> i+1 which together with x^ y i , 1^, x^ and 
y^ j are used to obtain the parameters for the ith way point. The heading of 
the straight-line segment between way points i - 1 and i is iJk , given by the 
relation 


• . y i " y i-l 
iK = arc tan 

1 i ' x i-i 


(-180 c 


< < 180 


( 2 ) 


The heading change ^ in the circular segment near way point i is 

(- 180 ' < * ti s 180 ') 

v ti TscT vi+i *i) 


(3) 
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where the direction of the turn is to the right for ^ > 0 and to the left 
for 4^ < 0. Next, calculate the quantity Ik shown in figure 3. 


b. = R. tan 
x i 


ti 


(4) 


The length of the straight segment d^ +1 is then 


d i.r - h 


(h.l * ®) 


(5) 


If d^ +i from equation (5) is less than zero, adjacent turns overlap, the 
computation of the trajectory must stop, and the pilot is given a diagnostic 
message, such as "way points i and i + 1 are too close." The coordinates 
for the end of the turn and the center of the turn are 


x Qi - x i * b i cos *i.l 


( 6 ) 


y Q1 - y 4 * bj sin * 1+J 


(7) 


x Ri = x Qi ; R i Sln h.l 


( 8 ) 


y Ri * y Qi * R i cos *1.1- (9) 

where the upper sign is chosen for a right turn and the lower sign for a left 
turn. Next, the distance l ^ is calculated: 


h c ^/( x i-i ' x i) 2 * ( y i-i - y i) 2 - b i (h 2 °) 


CIO) 
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If l ^ is less than zero, the calculation cannot continue and the pilot is 
given a diagnostic message such as "way points i - 1 and i are too close." 

If l j £ 0, the iteration is completed by calculating the coordinates of the 
beginning of the turn: 

Xpi ~ x i “ b i cos ** (H) 

y Pi = y i " ^i sin ^i O 2 ) 

These iterations are continued until the first way point is reached. Since 
the way points entered do not always yield a flyable trajectory, onboard 
calculation of the trajectory generally will require a system that permits the 
pilot to correct errors after the system has issued a diagnostic message. 

Final Heading Way Points 

Final heading way points are illustrated in figure 2 by way points 3, 6, 
and 7. Instead of rounding the corner at the intersection of two lines, the 
trajectory for this type passes through the way point at the instant the turn 
toward the next way point has been completed. Thus, the aircraft begins its 
flight along the straight-line segment exactly over the way point. There are 
two reasons for introducing this type of way point. First, it simplifies the 
specification of some trajectory segments, such as the turn at way point 3, 
which contains more than 180°. Recall that the turn at an ordinary way point 
is limited to less than 180°. Second, this type is required if the arrival 
time at the way point is specified. Specification of arrival time at an 
ordinary way point lacks precision since the way-point coordinates themselves 
do not fall on the trajectory. By requiring all aircraft that are merging to 
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fly through a point on the merging path with a common heading, the assignment 
of arrival times at the merging way point can be used to achieve precise 
spacing of aircraft. 

In the specification of a trajectory, ordinary and final heading way 
points can be alternated in arbitrary fashion. The general procedure for 
calculating the trajectory parameters is illustrated in figure 4 with a final 
heading way point i embedded between way points i - 1 and i + 1 of 
arbitrary type. As before, the trajectory is computed backwards from the 
last way point. The final heading to be achieved at way point i was 

previously determined in the calculation for way point i + 1 and is there- 
fore a known quantity. It is evident from figure 4 that the desired final 
heading at way point i can be achieved with two trajectories, one ending 
with a left turn, the other with a right turn. The criterion for selection 
is to choose the one with the shorter path length between way points i - 1 
and i. To make the selection, the coordinates of the centers of the two 
turns are calculated: 


*Ri ’ x i - R i si ” Hu 


y Ri ' y i * R i cos *i.l 


x Li ■ x i * R t sin *i-i 


y. . = y. - R. cos ip. 

7 Li 7 x x l+l 


(Right turn) 


(Left turn) 


(13) 

(14) 

(15) 

(16) 
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Then, the distances squared from way point i - 1 to each center are 



(17) 

(18) 


From the geometry of the construction in figure 4, it can be seen that the 
trajectory with the shorter path length also has associated with it the 
shorter of the two distances d R ^ and d^. Thus, the right-turn trajectory 
is chosen if d R ^ is greater than d^. If the way point i - 1 lies on 
the line determined by way points i and i + 1, both trajectories have the 
same length and the direction of the turn, if a turn is required, must be 
selected on the basis of another criterion. The next step is to determine 
whether way point i - 1 lies inside or outside the circle used to define 

the turn. Suppose the right turn was previously selected; then the way 

point lies inside if d„. < R? and outside or on the circle if d^. > R?. 

If this former case is true, the trajectory is not feasible, further 
computation of the trajectory stops, and the pilot is given the message 
"way points i and i - 1 are too close." In the latter case, the calcu- 
lation continues with the computation of the heading i^ R ^ of the directed 
line d R ^ if a right turn is required or the heading of the directed 

line d^ if a left turn is required: 


, y Ri y i-l 

’Vi = arc tan x ITT 

dKi x Ri x 1 . 1 


(Right turn) 


(19) 


*dLi = arC tan x 


y Li - y i-l 
Li " X i-1 


(Left turn) 


( 20 ) 
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Next, the angle is computed for a right turn: 


R. 

i/j R . = arc sin -yi- 
Rl d Ri 


( 21 ) 


or for a left turn: 


R. 

, • i 

i^ T . = arc sin - 3 — 

Ll d u 


( 22 ) 


The heading i/k is then 


1 mod . . 

^i * 180® ^dRi " ^Ri 


(23) 


if the trajectory contains a right turn and 


, mod , 
h = w *dLi - *U 


(24) 


if it contains a left turn. The length d^ of the straight-line segment from 
way point 1-1 to the beginning of the turn is 



(Right turn) (25) 


(Left turn) (26) 


Finally, the heading change ^ ^ in the turn is obtained from the difference 
between and 1 , and the coordinates of the beginning of the turn to 
way point i are 
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(27) 


x Pi - x i-l * d i cos *i 


y n . = y. , + d. sin ip. 
'Pi ' l-l l y i 


(28) 


Technique of Flying to the First Way Point 

In the preceding discussion, the trajectory from the first to the last 
way point was synthesized. To complete the synthesis, the trajectory from 
the aircraft's initial position and heading to the first way point must be 
constructed. 

The construction procedure depends on the type of the first way point. 

If it is an ordinary way point, the trajectory is constructed on the basis of 

the rule that the aircraft will turn from its current heading toward the first 
way point in the direction that minimizes the total path length to the first 
way point. This 'criterion is the same as that used in constructing the tra- 
jectory from way point i + 1 to way point i, where way point i is of the 
fixed final heading type. Thus, the procedures of the preceding section apply 
if the initial heading is identified with iJk + 1 , the initial position is 

identified with the coordinates of way point i, and the coordinates of the 
first way point are identified with those of way point i - 1, and the airplane 
traverses figure 4 in the opposite sense. 

If, however, the first way point is the fixed final heading type, a pro- 
cedure different from those discussed so far must be used. This guidance 

problem can be stated as follows: determine a trajectory that starts from a 

given initial position and heading and leads to a position with a specified 
final heading. Since this problem has been dealt with extensively in 


- 14 - 



references 1, 2, 3, 6, and 7, only a brief discussion of results is given 
here. References 1 and 2 give a minimum time solution to this problem, 
reference 2 deriving the optimum control law. Simplified solutions are dis- 
cussed in references 3 and 6. In the simplified treatment, the required 
trajectory consists of a turn, a straight flight, and then another turn. 

The parameters of the three segments are chosen to satisfy the initial and 
final conditions of the problem. 

Altitude Profile 

The calculation of the altitude profile is simple and requires only the 
determination of a flight-path angle y^ between way points i - 1 and i. 
Since the horizontal path length between way points is known, the flight-path 


angle y^ is given by 


arc tan 



(29) 


where s^ is the horizontal path length between way points i and i - 1, and 
h^ and are the altitudes specified at way points i and i - 1, 

respectively. In order to compute s^ from the segments of turns and straight 
lines it is necessary to define at what point on the horizontal trajectory the 
specified way-point altitude is to be achieved. The rule used is that the way- 
point altitude must be achieved exactly at the end of the turn for a given way 
point. This rule is used for both ordinary and final heading way points. The 
last step in the altitude profile computation is to check whether each y^ lies 
within the range of permissible flight -path angles for the aircraft. 

It is possible and perhaps desirable to define more complicated altitude 
profiles between way points. Profiles that minimize a performance function 
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such as the fuel consumed could be valuable. However, the limited size of 
the airborne computer to be used in flight tests of this guidance system, 
together with the requirement for onboard computation of the trajectory, does 
not permit consideration of more complex techniques at this time. 

SPEED PROFILE COMPUTATION 

Precise time control of aircraft is achieved by determining a feasible 
speed profile along the already computed 3-D trajectory. A speed profile is 
said to be feasible if it satisfies the following conditions: 

1. The airspeed remains between the minimum and maximum airspeed 
restrictions imposed along the trajectory. 

2. The rate of change of airspeed does not exceed the acceleration/ 
deceleration capabilities of the aircraft. 

3. The resulting ground speed yields the desired arrival times at those 
way points where' they have been prescribed by ATC (such way points are 
referred to as time-controlled way points). 

The first two conditions imply that the flying time between any two points on 
the trajectory is bounded above and below by the minimum and maximum times 
corresponding to the maximum and minimum airspeeds, respectively. Consequently, 
if arrival times at the time^controlled way points are assigned arbitrarily, 
then a feasible speed profile may not exist. In order to insure the existence 
of a feasible speed profile, ATC assigns arrival times based on the minimum 
and maximum possible flying times between successive pairs of time-controlled 
way points. 
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The determination of a feasible speed profile in the absence of wind 
would be a relatively simple task, for in this case airspeed and ground speed 
are identical. In the presence of wind, however, the relationship between 
airspeed and ground speed is highly nonlinear, and the problem becomes con- 
siderably more complex. By assuming a steady wind, the expression for the 

magnitude of the ground speed V (t) at any time t is given by 

S 

V g Ct) =^(t) - V* sin 2 * w Ct) + V w cos * w (t) (30) 

where V (t) and V are the magnitudes of the airspeed and windspeed, 
respectively, and i|» w (t) is the angle fromthe wind direction to the ground 

The differential equation governing 

(Straight flight) (31a) 

(Circular flight) (31b) 

where R is the radius of turn. Exact analytic expressions for ip (t) and 

V (t) can be found only in the case of straight flight with constant airspeed. 
S 

A considerable simplification can be made by assuming that for all t 

/ v \ 2 

w ] « 1 (32) 

\ v . (t v 

Supporting evidence indicates that inequality (32) is a good approximation 
not only for CTOL but STOL aircraft operations as well. Under this assumption 
equation (30) can be written as 

V g (t) = V a (t) + V w cos * W C t) (33) 


heading, measured positive clockwise. 
^ w (t) is: 

<W w (t) 
dt 


0 

I V l) 
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Before considering the speed profiles for straight and curved flight 
in more detail, it is necessary to establish certain desirable character- 


istics for airspeed profiles. In the design of the 4-D guidance system 
described in this paper, the point of view was adopted that the airspeed 
profile should be a piecewise linear function of time, that is, 

V a (t) = V a( t k) + a k (* " *k) 0* = * = H+l' k = °» 1, * *) (34) 

where a^ is the constant value of acceleration/deceleration in the inter- 
val Furthermore, changes in airspeed should occur only at a 

few places along the trajectory, preferably at those points where the mini- 
mum and maximum admissible airspeeds change. These requirements were 
dictated by considerations of passenger comfort, pilot workload, and 
simplicity of implementation. A typical airspeed profile possessing these 
characteristics is shown in figure 5. 

Straight Flight 

Let the desired airspeed along a straight flight segment of length d^ 

be given by equation (34). Then the analytic expressions for , V , and 

" g 

k+l are 

*W< 1 > = *w(h) 

V'> - v a(h) * 


(35) 


a k( l ‘ l k) * V w cos ♦«(*) 


(36) 
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(37) 


-V 


k+1 


= t. 


M +jEH ~ 5 a 


Curved Flight 

Let a curved flight segment consist of a circular turn of radians 

with turning radius > 0 for a right turn, i|/ tk < 0 for a left turn). 

If the desired airspeed along the circular arc is given by equation (34) , then 
no analytic expressions can be found for < P w , V g , and t^ +1 . Since numerical 
integration for determining the speed profile would be prohibitive, a different , 
approach is needed. It so happens that equation (31b) becomes integrable if 
the airspeed has the following form: 


V 


,<« ■ v Jh) ♦•*(*♦ (* - *0 ' * 5 **>) 


(38) 


In this case the analytic expressions for i/> w , V , and t^ +1 are 


-1 


T w' 


C 1 

tan" 1 ^- tan 


\+ _ S ...... 

^ ^ t 2 ) 


7 ^- tan 

C2 

2 

/ 

/ 2R* sgn * tk 

( 2v *n [ 

i 


/ v cos * (t) \ / t _ t \ + v 

: (t) = MM + \{ 1 + - v a( t k ) ) V J w 


cos V^ w (t) 


l (39) 


(40) 
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where Ci, C 2 , and t^ are defined by 
C ‘" V a( t k)* V „ 

2R k sgn(* tk ) | Jc /* w (t k ) * * tk ) 

Vl Cl | tan [cl tan \ 2 / 

(Note that t^ +1 = t^ + j if a^ = 0.) Although given by equation (38) 

is not a linear function of time, for values of ^ a (' t ] c ) and satisfying 

inequality (32), V (t) turns out to be very nearly linear. This condition 

£L 

is illustrated in figure 6, which shows the airspeed and ground speed along 
a section of the example trajectory of figure 2. 

The only remaining quantity still to be determined is the desired air- 
speed profile. Since the earliest and latest possible arrival times are 
achieved by flying the aircraft on the boundaries of the admissible speed 
ranges, the actual airspeed profile corresponding to an intermediate arrival 
time must lie between the speed boundaries. (See fig. 5.) The nonlinear 
relationship between arrival time and airspeed necessitates the use of an 
iterative procedure for the determination of the desired airspeed profile. 
Basically, the procedure adjusts the cruising speed level between each pair 
of time-controlled way points so that the prescribed arrival times are 
achieved. This adjustment is made by using expressions of the form of '■ 
equations (37) and (41) . : 
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A final remark concerning the wind is in order. In this paper only the 
case of a steady wind is considered. It is well-known, however, that both 
the magnitude and direction of the wind are functions of the altitude. If 
these functions were known, they could be easily incorporated in the speed 
profile computation. 

AN EXAMPLE OF 4-D TRAJECTORY SYNTHESIS 

The preceding sections described the techniques used to calculate the 
two major elements of a 4-D trajectory, the 3-D path and the airspeed profile 
along the path. These elements must now be assembled to produce the complete 
reference trajectory consisting of the reference states (position, altitude, 
and heading) and the reference controls (turning radius, airspeed, and flight- 
path angle) as a function of time from initial time to final time. 

The reference states and commands are calculated by a procedure which 
makes use of the chosen parameterization of the 4-D trajectory to minimize 
computer storage. The calculation is done in two steps. In the first step, 
the 3-D path and the airspeed acceleration time history along the path are 
used to construct a command table consisting of a sequence of control inputs 
arranged in chronological order. Since the reference controls are piecewise 
constant in time, the command table gives the values of the reference controls 
only at time instants where they change to new constant values. In the second 
step, the reference states between command times are computed analytically 
from the initial condition at the command time and the value of the controls 
during the command interval. Compared with the technique of storing the 
reference trajectory at a large number of time instants, this technique uses 
significantly less storage, an important consideration in implementing the 
technique on an airborne computer. 
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The example trajectory shown in figure 2 is used to illustrate the tech- 
nique of 4-D trajectory synthesis described in this paper. The pilot specified 
the trajectory to the system by entering the data given in table I. In the 
onboard system the way-point types are replaced by numerical codes. Both the 
initial position (lift-off) and the touchdown point are treated as final head- 
ing way points in synthesizing the trajectory. The initial heading and the 
runway heading associated with these way points are both 0° in this case. 

Final heading way point 6 and the touchdown point are time control points with 
arrival times of 300 and 350 sec, respectively. The wind is assumed to be 
from 0° at 25 ft/sec. The airspeed range specified for each way point in 
table I is valid from the end of the turn performed at that way point to the 
end of the turn performed at the next way point. This procedure requires 
choosing starting times of decelerations/accelerations such that the airspeed 
will fall in the next speed range at precisely the end of the turn. These 
fixed boundary conditions are met by synthesizing the airspeed profile back- 
ward from the last way point. 

The 4-D command sequence generated for this input is given in table II. 

The minimum and maximum arrival times at the two time control way points 
(WP6 and WP7) are given in the headnote of the table. There are 17 command 
times in this example. Columns 3 to 8 give the states at the command times 
and columns 9 to 11 the piecewise constant controls between command times. 

Space limitations prevent giving the equations for computing the reference 
states between command times. For the same reason the equation for computing 
the instantaneous bank angle from the airspeed, heading, wind vector, and the 
turning radius is not given. Bank angle and flight-path angle commands are 
applied to the aircraft stability augmentation system (SAS) slightly in advance 
of those given in the table to minimize errors due to the finite rates of these 
quantities. Computation time of this trajectory on the IBM 360 is 0.5 sec. 
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PERTURBATION EQUATIONS AND CONTROL LAW 

A synthesis procedure having been developed for the 4-D trajectory, which 
will now be referred to as the reference trajectory, the next step is to design 
a control law for flying it. The design of a control law for this problem, 
which is based on the technique described in reference 7, is accomplished by 
means of a perturbation method. Design of the control law for the altitude 
channel will not be considered here since this channel is simple with minimal 
coupling to the other channels. 

The nonlinear dynamical equations, from which the perturbation equations 
are derived, are as follows: 

X = V cos T : 

a (42) 

Y = V a sin V (43) 

• _ g(tan $) (44) 

* V 

a 

where V is the airspeed, X and Y are components of V , ! is the head- 

cl 3. 

♦ 

ing, 4' the heading rate, $ the bank angle, and g the acceleration of 
gravity. The wind is assumed to be zero and the flight -path angle y small 
so that cos y ~ 1 . 

The perturbation equations are obtained from the nonlinear equations by 
expansion in a Taylor series about a moving target reference system as illus- 
trated in figure 7. The origin and positive X-axis of this system at any given 
time are the reference position and the direction of the aircraft flying the 
reference trajectory, respectively. In figure 7, X r> Y , and 4* r refer to 
the reference position and heading and X fl , Y & , and 4^, to the aircraft 
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position and heading in the runway- centered coordinate system. The linear 
differential equations obtained from a Taylor series expansion in the per- 
turbed quantities X, Y, and Y shown in figure 7 are 


x = v + fy 1 
a r J 

y = v ip - Y x 
1 r r 


(45) 



(46) 


where v = V - V and ' <J> » $ - $ . The advantage of deriving the perturbation 

equations in this target referenced system is that terms involving sin and 
cos Y , which would otherwise appear, are eliminated; thus the perturbation 
equations for curved trajectories are simplified. Equations (45) and (46) also 
show that x and -y are coupled when j* 0, whereas the gain of the if< 
channel is inversely proportional to V r ; thus, Y r and are parameters 

which depend on the reference trajectory. A control law for nulling the per- 

' ♦ 
turbed quantities will therefore have to contain and as parameters. 

The control variables of the aircraft for tracking the reference aircraft 
are the bank angle $ and the speed V . A linear model of the combined auto- 

d- 3. 

pilot and aircraft dynamical response for a bank angle and velocity command 
system can be approximated by the following equations: 


♦ ■ ^ (-* * V ' Vc) (47) 

\ ■ r- (-’'a * k v v a - k v v c) < 48 ) 
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where <J> c and v £ are the command inputs and <j> and v & are the response. 

The parameters in these two equations, k^, k y , x^, and x y , were deduced by 
matching the step responses of these equations to those of a currently 
in-service four-engine jet aircraft with an autopilot and autothrottle. 

Their numerical values are 0.375 sec” 1 , 0.167 sec -1 , 1.04 sec, and 4.17 sec, 
respectively. The following control law is chosen for nulling the perturbed 
quantities x, y, and \ p: 

*c ’ - k *xV - V - Vr* t49) 

v c = - k vx x - k vyV < 50 > 

Note that equations (49) and (50) contain and V r as parameters. This 

parameterization .has been. found to be effective in achieving acceptable per- 
formance of the control law for the class of reference trajectories of interest 
here . 

The governing factors' for determining the numerical values of the five gains 
k vx » k y y., k^ x , k^, and are (a) the accuracy of the navigation data, (b) 

the allowable bank angle and throttle activity for passenger comfort, and (c) 
the accuracy of following the synthesized reference trajectory. A root-locus 
analysis of the closed loop system indicates that a good compromise between con- 
flicting requirements (b) and (c) is to use 0.0002 rad/ft for k^, 0.004 rad/rad 
(ft/sec) for k^, 0.0001 rad/ft (sec/rad) for k^ x , 0.04 (ft/sec)/ft for k yx> 
and 0.15 (ft/sec)/ft (sec/rad) for k y y. This combination of gain constants was 
obtained by trial and error by using a root-locus analysis of equations (45) to. 
(50). The roots corresponding to this set of gains yield reasonable frequency 
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and damping for all values of ¥ from zero to 6°/sec. A root-locus plot of 
equations (45) to (50) as a function of ¥ for the choice of gains given here 
can be found in reference 7. 

SIMULATION RESULTS 

Figure 8 shows the block diagram of the simulation used to evaluate the 
guidance system. The general flow of computations in the trajectory synthesis 
algorithm is indicated inside the block drawn with dashed lines. The final 
product of the synthesis computations is the command schedule. The reference 
states and controls generator uses the command schedule, clock time, and the 
measured wind vector to compute the reference states and controls for each 
control time interval. An interval of 0.1 second was used in the simulation. 

By use of the measured wind vector, ground speed and ground heading are con- 
verted to reference airspeed V r and airspeed heading 4^. 

The control loop used to fly the aircraft along the reference trajectory 

is shown in detail. The first step is to compute the perturbed quantities x, 

y, z, and 4>, which are obtained in the transformations given in the bottom 

block. These quantities are multiplied by the appropriate gains and are added 

as required to form the perturbation controls v , <j> , and z . They are then 

subtracted from the reference controls to form the autopilot and autothrottle 

inputs consisting of the command bank angle $ , the command airspeed V , and 

c c 

the command altitude rate Z . 

c 

A simplified model of the autopilot, autothrottle, and aircraft dynamics 
was developed especially for use in 4-D guidance and air traffic control simu- 
lation studies. A detailed flow chart for this model is given in reference 7. 
The model consists of a tenth-order dynamic system with hard limits on roll. 
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roll rate, airspeed, airspeed acceleration, flight -path angle, and flight -path- 
angle rate. The actual wind vector is also an input to the model. Output 
quantities are the actual aircraft states. From these quantities the naviga- 
tion system simulation obtains the measured aircraft states and the estimated 
wind vector. 

A complete analysis of simulation results cannot be given within the 
length of this paper. Only the response of the control law to track a STOL 
type reference trajectory in the horizontal plane consisting of a 360° circu- 
lar segment with radius of 1220 ft and an airspeed of 135 ft/sec will be shown. 
The reference trajectory, which has a duration of 56.7 sec, is generated with 
two final heading way points located on the circle as shown in figure 9. This 
reference trajectory is a severe test of the control law since to fly it 
requires a reference bank angle of 25°, almost equal to the bank angle limit 
of 30° used in the aircraft simulation; this angle leaves little bank margin 
for nulling out errors. Figure 9 also shows the trajectories of the simu- 
lated aircraft for two initial conditions and an error in the wind estimate. 
The, position of the reference aircraft and of the simulated aircraft is marked 
every 10 sec along the trajectories. 

Starting from the two initial conditions, the simulated aircraft locks 
onto the reference trajectory after 30 sec of flight even though a period of 
bank angle limiting occurs (not shown in fig. 9) while the control law nulls 
the errors. To evaluate the effect of wind-estimate errors on tracking 
accuracy, an 8.45 ft/sec (5 knots) constant wind error was introduced. Nor- 
mally, a wind estimator, which is part of the navigation system, would 
observe this error to a degree and refine its estimate, but in this case the 
estimator was disabled. The resulting tracking error is 180 ft at the end of 
the trajectory and indicates the importance of accurate wind estimates in 
precision aircraft control. 
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CONCLUSION 


The chief advantage of the approach to 4-D guidance described here is the 
ability to specify and compute complex trajectories in flight. This feature 
is a highly desirable one from the pilot's viewpoint. Another advantage is 
that the technique is not strongly dependent on the aircraft type, since the 
only aircraft parameters used in synthesizing the trajectories are performance 
limitations, which are treated as parameters. Furthermore, the guidance tech- 
nique can be integrated with a ground-based scheduling technique to form a 
complete air traffic control system. The precision of trajectory control and 
arrival time achieved with the system provides a solid base for reducing 
separation requirements and increasing landing rates in future air traffic 
control system. 
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TABLE I. - INPUT QUANTITIES REQUIRED FOR EXAMPLE TRAJECTORY 

[initial heading, 0°; runway heading, 0°; airspeed acc/dec, 1.5 ft/sec 2 ; 
wind speed/direction, 25 ft/sec/O 0 ; initial airspeed, 135 ft/sec] 


Way-point 

number 

Way-point 

type 

Way-point 

coordinates, 

ft 


Turn 
radius , 
ft 

Airspeed 
range , 
ft/sec 

Time, 

sec 



X 

y 

fl 




1 

Initial position 

1,000 

• \ 

0 

0 

4000 

135 to 135 

0 

2 

Ordinary 

5,000 

0 

600 

4000 

203 to 304 

— 

3 

Final heading 

15,000 

0 

2000 

4000 

203 to 304 

— 

4 

Ordinary 

3,000 

-3000 

1500 

4000 

135 to 203 

— 

5 

Ordinary 

-1,500 

4000 

1000 

3000 

110 to 135 

— 

6 

Final heading 

-4,500 

4000 

1000 

1500 

110 to 135 

300 

7 

Final heading 
(touchdown) 

0 

0 

0 

1500 

110 to 110 

350 












1 TABLE II. - 4-D COMMAND SEQUENCE FOR EXAMPLE TRAJECTORY 

fr . /T to WP6 = 252/324 sec; T . /T to WP7 = 295/376 sec 
L min max min max 




States 

Controls 

Command 

sequence 

t. 

■1 

■■ 

h. 

4 >> 

V > 
g 

V 

v 2 

R, 

Y, 

number 

sec 

Km 1 

ft 

ft 

deg 

ft/sec 

ft/sec 

ft/sec 

ft 

deg 

1 

0 

1,000 

0 

0 

0 

110 

135.0 

mm 

Straight 

6.2 

2 

18.2 

3,249.6 

0 

243.2 

0 

137.7 

162.3 


4000 right 

6.2 

* 3 

39.6 

6,187.8 

1285.9 

600 

47.3 

174.1 

191.1 


Straight 

3.6 

4 

64.1 

9,383.6 

4749.1 

898.3 

47.3 

210.9 

227.9 

0 

Straight 

3.6 

5 

76.0 

11,091.6 

6594.7 

1057.5 

47.3 

210.9 

227.9 

0 

4000 left 

3.6 

6 

0 

143.7 

15,000 

0 

2000 

-166.0 

252.1 

227.8 

0 

Straight 

-2.0 

7 

* • 

159.8 

• 

11,071.9 

• 

-979.4 

1860.2 

• 

-166.0 

» 

252.1 

• 

227.8 

I 

-1.5 

Straight 

• 

i 

-2.0 

# 

« 

• 

14 

• 

300.3 

• 

-45,000 

• 

• 

0 

1 • 
590.0 

9 

0 

♦ 

• 

91.1 

• 

• 

116.1 

♦ 

• 

1.5 

• 

• 

Straight 

• 

• 

-7.5 

15 

300.5 

-4,482.5 

0 

587.7 

0 

91.1 

116.4 

0 

Straight 

-7.5 

16 

345.4 

-378.2 

0 

49.6 

0 

91.1 

116.4 

• -1.5 

Straight 

! -7.. 5 

.17 

350.0 

1 

0 

0 

0 

0 

85.0 

110.0 

0 

Straight 

0 
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Fig. 1. Typical terminal area route structure 
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Fig. 6. Portion of airspeed and groundspeed time history for 
example trajectory, wind = 25 ft/sec 
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Fig. 9. Simulated flight with initial condition errors 

and wind estimate error 


